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Abstract 

In this article, we take the tensor currents to interpolate the spin-singlet heavy quarkonium 
states h c and hb, and calculate the masses and decay constants using both the Borel sum 
rules and moments sum rules, the predictions Mh a = 3.461q'j3 GeV (or 3.471q \^ GeV) and 
M hb = 9.821^4 GeV (or 8.87lo'.io GeV ) are consistent with the experimental data. 
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1 Introduction 

i-C ' In 2011, the BABAR collaboration observed evidences for the spin-singlet bottomonium state 

h b (lP) in the sequential decays T(3S) -)■ n°h b (lP), h b (lP) -t yr) b (lS) pQ. Later, the Belle collab- 
oration reported the hrst observation of the spin-singlet bottomonium states h b {lP) and h b (2P) 
with the significances of 5.5 a and 11.2 a respectively in the collisions e + e~ — > h b (nP)iT + 7r~ at 
energies near the T(5S) resonance, and determined the masses Mh b (ip) = (9898.3 ± l.ll^) MeV 
and M hb( 2P) = (10259.8 ±0.61^) MeV 0. On the other hand, the mass of the spin-singlet 
charmonium state h c (lP) was updated from time to time since its first observation in the pp col- 
lisions by the R704 collaboration [3] , the average value listed in the Review of Particle Physics is 
M M1P) = (3525.41 ± 0.16) MeV @]. 

The QCD sum rules is a powerful theoretical tool in studying the heavy quarkonium states 
[HI!], and the existing works focus on the S'-wave heavy quarkonium states J/ip, r) c , T, i] b , and the 
■ P-wave spin-triplet heavy quarkonium states Xcj, XbJ, J = 0, 1, 2, while the works on the P-wave 

spin-singlet heavy quarkonium states h c and h b are few [6l [7] , where the interpolating currents 
Q{x)d fjL ^Q{x) are used, Q = c,b [6]. The heavy quarkonium states play an important role both in 
studying the interplays between the perturbative and nonperturbative QCD and in understanding 
the heavy quark dynamics due to the absence of the light quark contaminations. The P-wave 
spin-singlet heavy quarkonium states Hq have the spin-parity-charge-conjunction J PC = l" 1 , the 
tensor currents Q(x)a fM1/ Q(x) without derivatives have the following properties under the parity 
| and charge-conjunction transforms, 

Q{x)a^Q{x) A Q{x)a^Q(x) , 



> 



- 1—1 

X 



Q(x)a^Q(x) A -Q(x)a^Q(x), (1) 

where x M = (t,x) and x M = (t, — x). We also expect to interpolate the heavy quarkonium states 
\iq with the axial-vector currents Q(x)j tJ,r yc ) Q(x), which have the following properties under the 
parity and charge-conjunction transforms, 

Qixh^isQix) — > -<3(£)7m75<2(£) , 

Q{x)rroQ{x) A g(x)y 75 Q(x) . (2) 

The axial- vector currents couple potentially to the axial- vector heavy quarkonium states Xci and 
Xbi, which have the quantum numbers J PC = 1 ++ . In this article, we study the h c and h b with the 
tensor currents using the QCD sum rules, and make predictions to be confronted with experimental 
data. 
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The article is arranged as follows: we derive the QCD sum rules for the masses and decay 
constants of the heavy quarkonium states Jiq in Sect. 2; in Sect. 3, we present the numerical results 
and discussions; and Sect. 4 is reserved for our conclusions. 



2 QCD sum rules for the heavy quarkonium states Hq 

In the following, we write down the two-point correlation functions Tl^apip) in the QCD sum 
rules, 

IW(p) = t J d 4 xe^ x (0\T[j^(x)jip(0)}\0), (3) 

Jl„{x) = Q(x)a^Q(x), (4) 

where J^ u (x) — J^ u (x), J 2 v (x), the two interpolating currents are related with each other through 
the relation a^j 5 = i e ^P aafs . 

We can insert a complete set of intermediate hadronic states with the same quantum numbers 
as the current operators J^ v {x) into the correlation functions H^apip) to obtain the hadronic 
representation [6] . After isolating the ground state contribution from the the heavy quarkonium 
states Hq , we get the following result, 

P 

hW/3(» = , f2 Q n {-QuolPvPP ~ 9v[)P^Pa + Q^fiPvPa + gvaP^Pft) H , 

M h Q - p 

= n (p) (-g^aPvPp ~ QvliPuPa + g^PvPa + 9va.Pv.Pf}) H , (5) 

where the decay constants fh Q are defined by 

(0|^(0)|/iq(p)) = ]h Q {£ i ,Pv~e v p li ) 1 

(0\JU°)\ h QiP)) = ifh Q e^Xre x p T , (6) 
and the are the polarization vectors of the heavy quarkonium states Hq. 

We carry out the Borel transforms (and the derivatives) with respect to the variable P 2 = —p 2 
to obtain the Borel sum rules (and the moments sum rules), and obtain the following results at 
the phenomenological side, 



1 



r(n) V dp2 J 

1 f" Imn(s) _^ P hQ _<o 



n(».0 = x^(-^)V 2 )I^=w q£) 

1 f So ImII(s) /, 
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* Am" S ( S + 4m 2 Q 0" +1 = (Af, 2 Q +4m^)" +1 ' 1 j 

where the so are the continuum threshold parameters. 

In the following, we briefly outline the operator product expansion for the correlation functions 
II^i/ Qi 8(p) in perturbative QCD. We contract the quark fields in the correlation functions Tlnuaffip) 
with Wick theorem firstly: 



H-nuap(p) = i / d 4 xe lp ' x Tr{a fJ , u j 5 Sij(x)j 5 a a /3Sj i (-x)} 
'. J d i xe %p " x Tr {cT liV S i j{x)cr a pSji(—x)} , 
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where 
Sij(x) 



(2tt)' 



(Tke 



-ik-x 



9sG^a a PQt + mQ ) + $ + m Q )a a P 8 lJ (g 2 s GG) 



mQ 



(k 2 
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m Q k 2 + m Q $ + gsDcG^Qt + m QXfW* + fW)Qt + m Q ) 



(k 2 - m 2 Y 



(9) 

and t n = 4^-, the X n is the Gell-Mann matrix, the i, j are color indexes, and the (g 2 GG) is the gluon 
condensate [6], then complete the integrals both in the coordinate space and in the momentum 
space, finally obtain the correlation functions n M „ Q ^(p) at the level of the quark-gluon degrees of 
freedom. 

Once analytical expressions of the QCD spectral densities are obtained, then we can take the 
quark-hadron duality and perform the Borel transforms (and the derivatives) with respect to the 
variable P 2 = —p 2 to obtain the Borel sum rules (and the moments sum rules): 
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(10) 

1 x 3 + (l-x) 3 
dx — 

x(l — X) 

lQa 2 s (qq) 2 

81 



dx 



(n + l)(n + 2) 



2x(l — x) ) (toq 
{n + l)(n + 2)(n + 3)m 2 Q 
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(11) 



where 



P(s) 



3 
2^2 



dxx(\ — x) 



,a s GG 



12T 4 



dx 



x 3 + (1 - xf 



x(l — x) 



8(s-m 2 Q ) 



1 ,a s GG 



12T 2 
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x(l-x)> f 



2x(l - x) 
1 - 



) / dx5(s — rhg) - 
[x 2 + (1 - x) 2 ]s 

and x; = i 
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5(s -m 2 Q ), 
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In this article, we take 



into account the gluon condensate ( asGG ) and the four quark condensates (qq) 2 , and neglect the 
three-gluon condensate due to its small value. In calculations, we have used the equation of motion, 
= Tlq = u d s 9s(jlfit a q, and taken the approximation (ss) = (qq). 
We can eliminate the decay constants fh Q and obtain the QCD sum rules for the masses of the 
heavy quarkonium states hq, 

ds d(-l/T^) P( s ) ex P , 

(12) 
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(13) 
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3 Numerical results and discussions 



From the experimental data M hb{1P) = (9898.3 ± l-l±i;?) MeV, M, lb(2P) = (10259.8 ± 0.6t^) MeV 
[2], M hc(1P) = (3525.41 ± 0.16) MeV, M Xc2(2P) = (3927.2 ± 2.6) MeV jt], we obtain the threshold 
parameters = (16 ± 1) GeV 2 and = (104 ± 2) GeV 2 approximately. The vacuum conden- 
sates are taken to be the standard values (qq) = -(0.24 ± 0.1 GeV) 3 and (^f^-) = (0.33 GeV) 4 
at the energy scale n = 1 GeV [6j |8] . The Q-quark masses appearing in the perturbative terms 
are usually taken to be the pole masses in the QCD sum rules, while the choice of the toq in 
the leading-order coefficients of the higher-dimensional terms is arbitrary [S]. The MS masses 
to q( to q) relate with the pole masses niQ through the relation 

m Q {ml) = rn Q {l + —-Ri + ...\ . (14) 

In Fig.l, we plot the /iq masses Mh q with the variations of the Borel parameters T 2 and moments 
parameters n both for the MS masses mQ(mg) and pole masses itlq. From the figure, we can see 
that taking the MS masses m c — m c (m 2 c ) — 1.2 GeV and m b = rn b (m 2 ) = 4.2 GeV underestimates 
the masses Mh Q - In this article, we take the pole masses m c = m c = (1.4 ± 0.1) GeV and 
m b = m b = (4.7 ± 0.1) GeV. If we take the Borel parameters as T 2 = (8 - 10) GeV 2 and (12 - 
15) GeV 2 in the channels h c and h b respectively, the pole contributions are about (47 — 55)% and 
(49 — 58)% respectively, it is reliable to extract the ground state masses. On the other hand, 
the main contributions come from the perturbative terms, the operator product expansion is well 
convergent. 

In Fig. 2, we plot the masses Mh Q with the variations of the parameters n and £ in the moments 
sum rules. From the figure, we can see that if we take the values £ > 12 and > 25 in the channels 
h c and h b respectively, the variations of the moments parameters n lead to small uncertainties, 
which can be neglected safely. In this article, we take the values £ = 14 and 30 in the channels 
h c and h b respectively, which correspond to the effective parameters £ e // = 2.3 and 5 respectively. 
The effective parameters £ e // are defined as £ e // = — x)) = |, the (x(l — x)) denotes the 

average value, i.e. 



1 [x(l - x)] k [x{l-x)] h 



m Q + 4m Q^(l -x) m 2 Q + 4m|Ce// 



(15) 



Small £ e f f cannot warrant the convergence of the operator product expansion, and therefore stable 
sum rules with respective to the variations of the moments parameters n. 

Taking into account the uncertainties of the input parameters, we obtain the values of the 
masses and decay constants of the heavy quarkonium states h c and h b , which are shown in Figs. 3- 
4, 

M hc = 3.46±°;} 2 GeV, 
M hb = 9.82t°;i 2 GeV, 
fh a = 0.38^GeV, 

fh b = 0.38±°;E£GeV, (16) 



from the Borel sum rules, and 



M ha 


= 3.47±°;} 3 GeV, 




M hb 


= 8.87^;i°GeV, 




he 


= 0.38±8;g| GeV, 




fh b 


= 0.39+°;°«GeV, 


(17) 
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Figure 1: The masses Mh Q with the variations of the Borel parameters T 2 and moments param- 
eters n, where £ = 14 and 30 in the moments sum rules for the h c and hb, respectively. 
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Figure 2: The masses Mh Q with the variations of the parameters n and £ in the moments sum 
rules. 
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Figure 3: The masses M^ Q with the variations of the Borel parameters T 2 and moments param- 
eters n. 
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Figure 4: The decay constants fh Q with the variations of the Borel parameters T 2 and moments 
parameters n. 
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from the moments sum rules, the masses from both QCD sum rules are consistent with the exper- 
imental data, while the values from the moments sum rules are better. 



4 Conclusion 

In this article, we take the tensor currents to interpolate the P-wave spin-singlet heavy quarkonium 
states /iq, and calculate the masses and decay constants using both the Borel sum rules and 
moments sum rules, the numerical values of the masses are consistent with the experimental data. 
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